Coulomb blockade in metallic grains at large conductance: beyond the 

zero-dimensional limit 



(N 

O 

o 

C 



(N 






cd 



13 

C 

o 

(N 
> 

oo 

o 

O 

o 



c3 



i 

c 
o 
o 



X 



I.S. Beloborodov and A.V. Andreev 

Bell Laboratories, Lucent Technologies, Murray Hill, NJ 07974 

Department of Physics, University of Colorado, CB 390, Boulder, CO 80390 

(February 1, 2008) 

We consider mesoscopic fluctuations of the Coulomb blockade oscillations in a disordered metallic 
grain in an external magnetic field coupled to a metallic lead by a tunneling contact with a large 
conductance gr and capacitively coupled to a gate. For a finite mean level spacing 8 in the grain and 
at temperatures T>i5we obtain the oscillatory part of the correlator of thermodynamic potentials 
at different values of the gate voltage. The correlations decay algebraically with the gate voltage 
difference in contrast to the exponential decay that would result from the mean-field treatment of 
the Coulomb interaction. The results are valid beyond the zero-dimensional limit and describe the 
crossover between the unitary and orthogonal ensembles. 



PACS numbers: 73.23Hk, 73.40 Gk, 73.21.La 

Coulomb blockade is one of the most striking manifes- 
tations of electron interactions at low temperatures [Q. 
It can be observed by measuring, say the charge of a 
small metal particle connected by a tunneling contact to 
a metallic lead and capacitively coupled to a metallic gate 
which is maintained at voltage V g , as in Fig. fn. 

In the weak tunneling regime the charge in the grain 
is nearly quantized at low temperatures and exhibits a 
characteristic step-like dependence on the gate voltage 
V g . This leads to the Coulomb blockade oscillations in 
the gate voltage dependence of various physical quanti- 
ties. The amplitude of these oscillations decreases with 
increasing transparency of the contact and depends not 
only on the total conductance of the contact but on the 
individual transparencies of the tunneling channels. 

In the strong tunneling regime most of the studies of 
the Coulomb blockade concentrated on the limit of large 
grains with the vanishing mean level spacing S. In this 
limit the Coulomb blockade oscillations of the thermo- 
dynamic quantities vanish at perfect transmission P,p). 
In the case of a multi-channel tunneling contact the 
Coulomb blockade oscillations at large conductance have 
been shown to be exponentially small in the dimension- 
less tunneling conductance gx @||- The case of a dif- 
fusive random contact was considered in Refs. MM- At 
5 = the mesoscopic fluctuations of the Coulomb block- 
ade oscillations due to disorder in the grain vanish. 

The case of strong tunneling and finite S was studied 
in Ref. M where it was shown that the Coulomb block- 
ade oscillations remain finite even at perfect transmis- 
sion and exhibit strong mesoscopic fluctuations due to 
the presence of disorder in the grain. This study utilized 
the bosonization method of Ref. || which is applicable 
to contacts with weakly reflecting channels. 

Recently, much interest was focussed on granulated 
disordered metals PJHJ, which can be viewed as ex- 
tended networks of metallic grains in the strong tunneling 



regime. The low temperature regime, where mesoscopic 
fluctuations at finite S become observable is well within 
experimental reach JlQ[ |. It is therefore desirable to gen- 
eralize the study of charging effects at finite S and strong 
tunneling m to granulated metals. The formalism of 
Refs. p,B| is well suited to treat individual grains but is 
difficult to generalize to extended disordered systems. On 
the other hand such systems are conveniently described 
within the cr-model formalism. It is therefore important 
to incorporate the non-perturbative Coulomb blockade 
effects into the framework of the cr-model approach. 




FIG. 1. Schematic drawing of a disordered metallic grain 
coupled by a tunneling contact to a metallic lead. The charge 
of the grain is controlled by the gate voltage V g ~ g 



In this Letter we show that mesoscopic fluctuations 
of the Coulomb blockade oscillations in a metallic grain 
with finite 6 in the strong tunneling regime can be stud- 
ied using the cr-model. We concentrate on the thermo- 
dynamic quantities of a grain of size L with the Thou- 
less energy Et — D/L 2 , where D is the diffusion con- 
stant, placed in an external magnetic field H. The grain 
is connected to a metallic lead by a non-random multi- 
channel tunneling contact with the dimensionless con- 
ductance gr = 2irh/e 2 R 3> 1 where R is the resistance 
of the contact. It is also capacitively coupled to a gate, 
as in Fig. pi. We assume that the electron-electron inter- 



action in the lead may be neglected due to screening and 
treat the Coulomb interaction of electrons in the grain 
within the framework of the constant interaction model, 



H ( 



N 



(1) 



Here Ec is the charging energy, TV is the operator of the 
number of electrons in the grain, and q ex V g is the num- 
ber of electrons that minimizes the electrostatic energy 
of the grain. We consider the case Et < Ec, and de- 
scribe the crossover between the unitary and orthogonal 
ensembles. The case Et > Ec for the unitary ensemble 
was recently studied in Ref. |ll| by a different method. 

The single particle levels in the grain are broadened 
due to tunneling into the lead. We assume that tunnel- 
ing contact is broad, so that each state x in the dot is 
coupled to many lead states a. In this case the probabil- 
ity distribution P(T X ) of level half- widths T x is sharply 
peaked about the mean value T. We can therefore ne- 
glect the fluctuations of P(T X ) and consider them to 
be equal to the mean value T which can be expressed 
through the dimensionless conductance of the contact qt 
as r = £z_. The temperature is assumed to satisfy the 
conditions S < T < E C - 

The main results of this paper are the following ex- 
pressions for the oscillatory parts of the ensemble aver- 
aged thermodynamic potential (£l(q) ) OS c and of the ir- 
reducible correlator ((Q(q)Q(q'))) osc = (Q(q)Q(q r )) osc — 
(n(q)) osc (n(q')) osc &t5\q-q J \^>T + T, 



(n(q)) osc = ~E c \n 



Ec 
T + T 



cos(2irq), (2a) 



((n(qMq'))) osc = l -E 2 c K 



p=0,l 



? 2ni(q-q') g 



(2b) 



Here Ec = 9t i° e~ gT / 2 and B p can be conveniently ex- 
pressed through the quantity a p = iS(q — q')/2 + p/th, 
with the cooperon gap 1/th = DeH/hc, 
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(2c) 



Here the index p labels the diffuson, p = 0, and the 
cooperon, p = 1, degrees of freedom and d — 2, 3 is the 
dimensionality of the grain. 

The results (0) were obtained with logarithmic accu- 
racy in the one instanton approximation and are valid for 
T > g T E c hx(E c /T) exp(-g T /2). The Zeeman splitting 
of energy levels was neglected. Equation (Eah coincides 
with the result for the zero-dimensional grain pT| . At 
q — q'^f oo the function B p in Eq. (Eq) tends to zero, and 
the correlator (VL(q)Vt(q' )) osc factorizes; (VL(q)Vt(q' )) osc = 

(n(q)) osc (n( q ')) osc . 

Below we outline the derivation of Eq. (H). We con- 
sider an ensemble of disordered metallic grains in which 



electrons move in the presence of a random impurity po- 
tential V(r) which is taken to be white noise random po- 
tential with the variance (V(r)V(r')) = <5(r — r')/2irvr, 
where v is the density of states. The correlator of the 
thermodynamic potentials at different values of the gate 
voltage q and q' can be written as 



d 2 



(n(q)n(q>)) =T 2 -^—(Z"(q)Z a (q>)) \ 



<9a9a' 



(3) 



Here (Z a (q)Z a (q 1 )) is the ensemble averaged product 
of the replicated partition functions at different values 
of the gate voltage. For 8 <C T <C E c it can be ex- 
pressed as a functional integral over the Q-matrices, aris- 
ing from disorder averaging, and over the phases <pj which 
arise from decoupling the Coulomb interaction (gj) via the 
Hubbard-Stratonovich transformation, as follows: 

(Z a (q)Z a '(q')) = Y / (Z a Z a ') w e *-* , (4a) 

{W} 

„2 e ±2j 



(Z a Z a 



F[Q,cf\ = Tr 
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(4b) 
(4c) 



Equations (Q) were obtained by taking the integral over 
the static component of the field <j> by the saddle point 
approximation which is valid for T ^> 5 |12|. The re- 



sulting expression, Eq. (Eaf) , is written as a sum over the 
sets of winding numbers {W} m& which label distinct 
topological classes of the phase, <j>j{(3) = <f>j(0) + 2ttWj. 
The treatment of the cooperon degrees of freedom re- 
quires doubling the number of fcrmion fields. The co- 
variant derivative d in Eq. (J4q) is expressed through 
the time-reversal symmetry breaking matrix T3 ]14| as 
d = V + ieAr^/c, where A is the vector potential, and 
Tr denotes the trace taken over the spin, <j, coordinate, 
r, replica, i,j, and Matsubara, n, m, as well as the time- 
reversal indices. The matrix J in Eq. (Eq) is defined as 



,/ 



iqjS 



+ re-* (r) A TT /e* (T ' ) 



OijOfj (j' 



(5) 



where i = id T is diagonal in the Matsubara frequencies, 



£ = S nm e n = 5 nm TrT(2n+l), and A Tj . 



-ij sin[7rT(r- 



t')]. Finally, qj = q for 1 < j < a and qj = q' for 
a < j < a + a' . 

At frequencies smaller than the inverse mean free time 
1/r the effective low-energy theory may be written as a 
functional integral over the saddle point manifold of Q- 
matrices ]15[ . The low-energy degrees of freedom are usu- 
ally parameterized as Q — UQqU^ 1 , where the matrix 
Qo is diagonal in the Matsubara space, and represents 



the saddle point of the free energy F[Q,<f>] in Eq. ([k]) at 
<\> = 0. At frequencies below the Fermi energy it is given 
by Qo = A. For our present purposes, however, it is con- 
venient to write the Q-matrix as Q = TAT -1 , where A is 
the saddle point with the lowest free energy F[Q,(p] for 
a given <j> p3]. It is given by the equation 

d d P (.* {p- eJ ^Y *A X 



A = — 



(27T)° 



iJ + 



2m 



m + Yt 



(0) 



Note that since the operator J is coordinate-independent, 
the saddle point matrix A is uniform in space. This is a 
consequence of the constant interaction model, Eq. (||) . 

At 4>j(t) — 2itTtWj the matrix J is diagonal in the 
basis of Matsubara frequencies. Since it is Hermitian 
it may be diagonalized in some different basis for other 
configurations of 4>j(t) as well: J = S nrn 8ij8a- ta >J^, where 
n, m can be identified with the Matsubara indices only 
at 4>j(t) = 2-kTtWj. The saddle point A, Eq. (§) is also 
diagonal in this basis and in the low-frequency sector is 
given by A = sgne n S nm 5ijS atCr i , where e„ should be un- 
derstood as the Matsubara frequency corresponding to 
Jl at <£j(t) = 2ttTtW 1 . 

At T 3> 5 the integration over the soft modes of Q in 
Eq. ( pb| ) may be carried out by expressing the rotation 
matrices T in terms of the generators, T = exp(iP), and 
expanding the free energy F[Q,(j)} to second order in P: 



exp(i^) 



u — z 
1 — uz* 



u = cxp(27rzTr), 



(9) 



whereas for Wj — I' 4>j is given by Eq. (Bh with z — ► z' . 
The free energy F[X, <f>] = gx on the instanton configu- 
ration (0) and is independent of z and z' which param- 
eterize the widths and the positions of the instantons 
in replicas with Wj = I, I' and satisfy the inequalities: 
\z\,\z'\ < 1 for 1,1' = 1, and \z\,\z'\ > 1 for 1,1' = -1. 
We need to integrate over the fluctuations of cf> around 
the instanton configuration Eq. (m. To this end we write 
<j> = 4> z + <j>, where <j> represents the massive modes. The 
integration over the zero modes should be performed with 
the measure d 2 z/(l — \z\ 2 ) for \z\ < 1. For \z\ > 1 the 
inversion transformation z — ► 1/z again maps the inte- 
gration domain to \z\ < 1 with the integration measure 
d 2 z/(l — \z\ 2 ). All other fluctuations <\> have a large mass 
of order gx and can be integrated out in the gaussian 
approximation. Moreover, in carrying out this integra- 
tion one may replace the determinants A[</>] in Eq. (S) 
by their values A[^ z ] on the instanton configurations (9) 
due to their weak dependence on <f>. We then obtain 



1 + aa' \ e 



{Z a (q)Z« (gQ) 
(Z°(q)Z°'{q')) {Q} ,,,_,, 

d 2 zd 2 z'f{z)f{z') A {M>} [0 2 ] 
(l-\z\ 2 )(l-\z'\ 2 ) A {0} [0] ' 

|z|,|a'|<l 



2TTi(lq+l'q) 



(10) 



F[Q,ct>]=F[\, 
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p^0,l (T,<t' i : j;n,m 
- A - JL + P/th) 
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(7) 



The prime over the sum in Eq. (M) and below indicates 
that the summation goes only over such pairs n, m that 
e n > and e m < 0. The Gaussian integration over P in 
Eq. (flb|) then results in the following expression 



-{w} 



n n 

i,j;p p,m,n 



(4 - JL + Dp 2 + P /th) 



(8) 



where the power 4 arises from taking the product over 
the spin indices. 

The largest term in the sum over {W} in Eq. (Ea) has 
a monotonic dependence on the gate voltages, q and q' , 
and corresponds to {W} = {0} with all Wj = 0. The 
leading oscillatory contributions to this sum arise from 
the terms with {W} = {1,1'}, with 1,1' = ±1, having 
only one non-zero Wj — I for 1 < j < a and only one 
non-zero Wj = I' for a < j < a + a'. The dissipative free 
energy F[X, (j)] in Eq. (m) has a minimum when <j)j = 
in the replicas with Wj = 0. In the replica with Wj = I 
the instanton phase configuration [ [17[ which minimizes 
the free energy may be written in terms of a complex 
variable z M, 



where f(z) is given by the following expression JTT| 

2 i^ f 2tt 2 T\z\ 2 \ 

gT- E cexp|- gc(1 _' |z l |2) j 



/(*) 



2tt 3 T [1 



ttT 



\z\* 



(11) 



With the aid of Eq. (0) the determinant A^; .z/}[</> 2 ] in 
Eq. ( ^0| ) can be expressed through the eigenvalues of the 
operator J on the instanton configuration (ph. To this 
end it is convenient to express the operator J, Eq. (jq), 
in terms of the complex variable u in Eq. (H) , 



it/d 
~2~ 



2ttT 



= -ud„ 



iT 

2^ 



du' 1 



u u - 



u'l 



u z* 
(12) 



The integral in this equation is taken over the unit cir- 
cle, |it'| = 1, and the operator J is acting in the space of 
functions spanned by the fermionic Matsubara frequen- 



cies u 



-n-l/2 



In the case of the instantons, 1,1' = 1, 
(anti-instantons, 1,1' = —1) the functions u~™ -1 / 2 with 
negative (positive) Matsubara frequencies, are eigenfunc- 
tions of J with eigenvalues J n = e n — T — iqS/2. All other 
eigenvalues J n can be found from the following equa- 
tion ill 



1- 



E 

k=0 



2T\z 



2 k- 



T(Jn + iqS/2)+e k + T' 



(13) 



where the top (bottom) sign in the right hand side de- 
scribes the instanton (anti-instanton) case. For 5\q — 
q'\ ^> T + T the integral in Eq. ( J10| ) is dominated by 



and 



the region of short instantons, 1 - \z\ 2 <C ^2_ 
1 — |z'| 2 <C j^p- In this case the solutions of Eq. (fi~3) 
are given by J n = -iSq/2 ± [e n + T - 21>| 2n (l - \z\% 
n > 0. Next, we use these eigenvalues to evaluate the de- 
terminant in Eq. (ph. Since the second term in Eq. (IIQ) 
is multiplied by aa 1 , to obtain the correlator in Eq. (||) 
we only need to evaluate the ratio of the determinants 
in the second term in the limit a, a' — *■ 0. It is not dif- 
ficult to show, see appendix C of Ref. [pi, that in the 
limit a, a' -» 0, A {0} [(f> z ], A {M} [</> Z ] and A^!^}^] are 



equal to unity, and A/j j ; 

E E 



is given by 



A 



isr^|^"(i- 



{1,-1} 



(En+em + Dp^+ap)^ 



— p P,P 



A 



{0} 



(14) 



{0} 



with a p defined below Eq. (2b). The ratio — - 

is given by the complex conjugate of Eq. (111). Using 
Eqs. Q, © and (§) we obtain Eq. (§). 

In conclusion we note that: i) Equation (0) was ob- 
tained with logarithmic accuracy, and is valid for tem- 
peratures 5<T« E C - To derive it we also neglected 
the instanton-instanton interactions, which is valid for 

T» 5 ^ c lnf^exp(- 5T /2)- 

ii) The results for the orthogonal and the unitary en- 
sembles are obtained from Eq. (||) by taking the lim- 
its th —> oo and th — > respectively. If the Zeeman 
splitting is taken into account the degeneracy of diffu- 
sons/cooperons in Eq. (Eq) corresponding to different 
spin projections of the particles is lifted: two of them 
remain unchanged, and the other two acquire a correc- 
tion i8(q — q 1 ) — * iS(q — q 1 ) ± i2fiH, where /i is the Bohr 
magneton. 

iii) The correlator of differential capacitances, or ther- 
modynamic densities of states, may be obtained from 
Eq. (Eq) by differentiating it with respect to q and 



■ln z 



q . the correlations tall on as \q — q \ = in \q — q \ 
at |q — q'\ ^ Et/S. Equations ( pq ) and ( pc[ ) demon- 
strate the inadequacy of the mean-field treatment of the 
Coulomb interaction for the description of the Coulomb 
blockade oscillations. Under such treatment the Coulomb 
interaction in the grain would merely result in the lin- 
ear g-dependence of the self-consistent chemical poten- 
tial for the non-interacting problem. In contrast with 
Eq. (0) the corresponding oscillatory part of the cor- 
relator of densities of states would decay exponentially 



oc |g- q'\ 4 exp< -■',; 
numerical coefficient 



(%?!) 



d/2 



L8j, where 7^ is a 



iv) The mesoscopic fluctuations of the Coulomb block- 
ade oscillations in the strong tunneling regime were first 
studied in Ref. ^ . The single particle properties of the 
dot in this approach are described in terms of the scat- 
tering matrix at the contact. The long range correlations 
obtained in Ref. ||] resulted from neglecting the depen- 
dence of the scattering matrix on the gate voltage. 

v) In the CT-model approach the Coulomb blockade 
oscillations are described by the instantons of the in- 
field J17f. Their mesoscopic fluctuations are described 
by the fluctuations of the Q-matrix about the instanton 
configurations. The problem therefore reduces the eval- 
uation of diffusons on the instanton configurations of the 
0-field. 
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